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Abstract. We express the zeta function associated to the Laplacian operator on 
S 1 ' x M in terms of the zeta function associated to the Laplacian on M, where M 
is a compact connected Riemannian manifold. This gives formulas for the partition 
function of the associated physical model at low and high temperature for any compact 
domain M. Furthermore, we provide an exact formula for the zeta function at any 
value of r when M is a D-dimensional box or a Z?-dimensional torus; this allows a 
rigorous calculation of the zeta invariants and the analysis of the main thermodynamic 
functions associated to the physical models at finite temperature. 
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1. Introduction 



The zeta function regularization procedure is one of the most elegant and rigorous 
methods to deal with path integrals in quantum field theory. A great effort has been done 
in the last years to produce exact calculations in different cases of interest PH] [H] 
[T2] [H] with particular attention to exact calculation of the heat kernel coefficients [3]. 
Notice also some rigorous calculations appeared in the recent mathematical literature: 
[Zj P (20j |2H] [2E] [IE]. The purpose of this note is to show the result of a rigorous 
application of the zeta function regularization technique. In particular, we analyze the 
basic example considered by Hawking in [THj : a scalar field in a compact domain of 
the product space time S^ 2w x M at temperature T = -s, and provide a complete 
treatment of it. Our main motivation is that we can perform all the calculations in 
a rigorous way and determinate all the quantities of interest without introducing any 
approximation. Our main result is an effective formula expressing the partition function 
at finite temperature in terms of the geometry of the spatial background. This allows on 
one side to describe the behavior of the physical model at low and high temperature for 
any spatial domain (cfr with [S], where a similar analysis was performed for M = H 3 /T); 
on the other side, to provide explicit formulas that describe the model for some particular 
geometries, namely the D-dimensional box and the D-dimensional torus, at any value 
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of the temperature. Our technique is likely to be generalized to other situations, and 
there are works in progress in various directions. 

The partition function for a scalar field with action J[0] = — J <f)(x)A<j)(x)d(vol(x)) 
in the Euclidean space time and where A is a real elliptic self adjoint second order 
differential operator, can be formally described by the Feynman path integral Z = 
J £>0e p/M = [det(pA))~2 where p is some renormalization constant with dimension of 
mass or inverse length introduced by hand in order to obtain a physically consistent 
model [16jf. With periodic boundary condition in the imaginary time with period 
(3, the background geometry of the flat space time is described by the product space 
^0/2n x with product metric l(Bg, where (M,g) is a complete Riemannian manifold 
of dimension D; in particular, we assume M to be compact connected. Following 
Hawking j!6j . the path integral can be given a rigorous interpretation in terms of some 
zeta invariants associated to the underlying geometry. We introduce the zeta function 
associated to the operator A, ((s; A) = J2 x <=Sp A ^~ S > (where Sp A denotes the non 
vanishing part of the spectrum of A); then, defining the regularized determinant of A 
m H by logdetA = - £C(s;A)\ we get \ogZ = §C'(0;A) - |logpC(0; A). Under 
the identification T = 4, the partition function for the quantum theoretical model 
corresponds to the one for a canonical ensemble at temperature T. We will assume 
this point of view, and we will work out the partition function for a statistical 
system calculating the zeta invariants of the underlying geometry. This will allow us 
to introduce and analyze other interesting thermodynamic functions and to get useful 
information on the system described at finite temperature T. The main feature of this 
approach is that we can describe completely the partition function for any theory in 
the product space time at any value of the temperature, in terms of a zeta function 
at null temperature (that we will call geometric zeta function) depending only on the 
geometry of the space, namely of the background physical domain. To get this result, 
after introducing the zeta function associated to the model in the proper Euclidean 
setting, we use the particular form of the spectrum to decompose the zeta function. We 
easy get a function of the temperature T, smooth for all positive T, but the hard point 
is to get the right analytic extension at s — 0, in order to obtain the partition function. 
In particular, the analytic continuation is not uniform in T for small T, and hence the 
calculation of C'(0) m ust be performed at positive T. We are able to get this result, 
and we can do it uniformly in T for T in any closed interval of the real positive axis 
(Proposition El of Section EJ). Moreover, the result is effective, and we can deduce from 
it both the behaviors for low and high temperature. 

Studying the models at low temperature, we show that periodic boundary conditions on 
the spatial domain give an anomalous behavior for the zero rest mass scalar field that 
can be corrected by adding a non zero mass term, while studying the behavior of the 
pressure of the radiation as a function of the volume at different fixed temperatures, 
we prove the existence of a minimal non vanishing value of the volume below which 

| Beside in this work we will not treat renormalization problems, we will give effective formulas for Z 
for any p. 
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the force becomes attractive, as expected because of the Casimir effect |2| [22j- More 
precisely, the expansions at low temperature of the thermodynamic functions show the 
following effect of the boundary conditions on the physics of the model: the presence 
of a spatial zero mode produces an anomalous behavior of the entropy and the specific 
heat at low temperature. In particular, on a closed (compact, connected, with no 
boundary) geometry, we always have a zero mode for the Laplacian, due to the constant 
eigenfunction, and this produces a log T term in the logarithm of the partition function, 
and hence a logarithmic divergence in the entropy and a non zero specific heat at low T. 
Such a term disappears when a boundary is present and opportune boundary conditions 
are assumed. Furthermore, the presence of a non zero mass term cancels the logarithmic 
term and recasts a correct behavior for small T on both the domains, independently from 
the boundary conditions. This is in agreement with the appearance of a logarithmic 
divergence (infrared divergence) for a massless field classically (see for example |23j). 
Particularly meaningful in this context is the case D — 1, whose associated zeta function 
appears as the zeta function for the Laplacian with a constant potential on a cylinder 
or on a torus. Its s-expansion near s = is well known since Eisenstein and Kronecker 
|H2] , and relates to the theory of modular functions and L-series. The coefficient of the 
linear term can be expressed using the Dedekind eta function, and using the modular 
property of the last we can easily relate the behavior at low and high temperature of 
the partition function of the physical model. It is clear why the presence of a mass term 
cancels the singular behavior at low T: in fact, such a non homogeneous term breaks 
the modular property of the eta function. This will be discussed at the end of part 14.31 
The D dimensional version of the eta function appearing in the corollary of Proposition 
IHJ has not been studied yet. It does not exhibit any periodicity, but possible modular 
properties should be very important. A deeper analysis of such function should be very 
interesting both under a theoretical and applicative point of view. 
The study of the thermodynamic function at finite temperature allows to analyze other 
interesting phenomena. In particular, we investigate the pressure of the radiation at 
finite temperature as a function of the volume, and we show the existence, for any fixed 
temperature, of a critical volume Vq where the pressure of the radiation changes sign, 
becoming attractive. The analysis performed in part 14.51 of Section H] shows that below 
Vo the pressure is attractive and decreases like V^^d 1 , as expected by the Casimir effect 
[2| [22]- We also show that, up to renormalization, the critical volume depends on the 

temperature by the law TqV q d = const, and the boundary does not affect this effect. 
The work is organized as follows. In Section |21 we briefly recall some basic information 
about the geometry of the Laplace operator on compact manifolds. This will lead us 
to define the zeta function that describes the physical models we want to study and 
to state our main result in Proposition |3] and its corollary. In Section El we prove 
asymptotic expansions at high and low temperature. In Section EJ we study the zeta 
function introduced in Section |21 for some particular geometries at finite temperature, 
and we provide formulas for the main thermodynamic functions. 
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2. The zeta regularized partition function 

We recall some known facts about the geometry of the Laplace operator on compact 
manifolds. Let (M, dM, g) be a compact connected Riemannian manifold of dimension 
D with (possibly empty) boundary dM. Let (Am,BC) be the Laplace operator built 
using the metric g and acting on the opportune space of functions with boundary 
conditions BC on the boundary of M. If the boundary of M is not empty, Dirichlet or 
Neumann boundary conditions can be chosen. Let Sp(— Am, BC) denotes the spectrum 
of the associated boundary value problem. Then, Sp(— Am , BC) is a discrete set of non 
negative (real) numbers: {\k}keK, K C N, and ker(— Am) is a finite set of rank fC. Let 
Sp denotes the set of the positive eigenvalues and K the set of their indices. If dM 
is empty, then /C = #ker(— Am) = 1, the constant function. In this situation, we have 

a 2 2 

the Weyl formula, that gives the behavior of the large eigenvalues, ~ - -rh~D , 

(volMvolB D ) D 

where B D is the unit ball in R D , and voI^d = [H] EH C3, and the heat 

kernel expansion ^2 k&K e~ Xkt = t" J2JLo e j^, when t — > + [12]. The formulas above 
hold also if a regular potential term is added to the Laplacian and eo = mXM D 1151 . 

(47r)"2" 

Before introducing the zeta function, we need the following technical lemma, that can 
be easily deduce from the Young's inequality. 

Lemma 1 For all real positive a and b, integer positive n and k, n a + k b > {nk)°^. 

We are now able to introduce the associated zeta function. We want to do this in 
slightly more general setting, namely we assume a possible mass term, and we consider 
the operator —Am + q, with real q > 0. The kernel of this operator depends on q, but 
its rank is discontinuous at q — 0. We will see in the following, that all the results 
we are able to prove, are true uniformly in q, for all q in some fixed closed interval of 
the positive real axis. This suggests to deal independently with the zero mass case. 
On the other side, to simplify notation and avoid to give always two formulas, it is 
preferable to choose the following alternative approach. Let's introduce the function 

fC„ = < ^ '^^^ , and assume q to be fixed. Then, we can write a unique formula 
I K if q — 0, 

for both the cases, q = and q ^ 0, that is clearly non smooth in q. Also, we state now 
once for ever that all the following formulas are smooth and uniform in q for q in any 
fixed closed interval of the positive real axis. 
We first introduce the zeta function at null temperature, 

C(s; 0, q) = C(s; -Am + q, BC) = ^'(A fc + q)'*, 

for Re(s) > y , where the notation means that the (possible) zero terms must be omitted 
in the sum; we will call this function the geometric zeta function associated to the model. 
Next, consider the product manifold S 1 x M, with the product metric d 2 t(Bg, Laplacian 
A = J^- + Am, and periodic boundary conditions u(0) = u{(3) on S 1 parameterized 
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by t G [0,(3]. Then, Sp(— A, BC) = {(2iiTn) 2 + \k}nez,keK, an d the associated zeta 
function is§ 

C(a; 2vrT, q) = ^' [(2vrTn) 2 + A fe + g]-, 

for Re(s) > by Lemma ^ The zeta regularized partition functions for a (possibly) 
massive scalar field whose underlying geometry is (M, dM, g) is 

logZ(T,g) = ~\ogdet(p(-A + q),BC) = ^'(0; 2ttT, q) - I logpC(0; 2vrT, q). 

As stated in the introduction, our approach is to study the zeta functions as 
mathematical objects, providing all the zeta invariants, and hence to write down 
formulas for the partition function of the physical models and to calculate the 
thermodynamic functions. In the most general setting of this section, all basic 
information about the zeta function can be obtained using classical methods (see ^Hj or 
We summarize them in Propositions ^ and |21 We also provide the fundamental 
analytic representation for the zeta function, Proposition EJ that will be the starting 
point for all the proofs in the following. From now on, we will use the variable y = 2nT 
to simplify notation. 

Proposition 1 The function £(s; 0, g) has an analytic continuation to the whole 
complex s-plane up to a set of simple poles at the value of s = for I = 

0, 1,2,.. that are not non-positive integers, with residua: Resi (C(s; 0, q), s — = 
r(i=i) ^j,k>o,j+2k=i ( "~W~ e 3 ( l k ' at zero and negative integers, -m = 0,-1,-2,-3,...: 
C(-m; 0, q) = E,- fc > 0J+2fc=2m+i > ^^g* - ^ . 

Proposition 2 The function ({s;y,q) has an analytic continuation to the whole 
complex s-plane, uniformly in y fory in any closed interval of the positive real axis, up to 
a set of simple poles at the value of s = , fori = 0, 1, 2, . . ., that are not non-positive 
integers, with residua\\: Resi y, q), s = ^ ± ^ d ) = r (£^i) \ ^j,k>o,j+2k=i ^S~ e 3l k > 
s = is a regular point and £(0; y,q) = & J2j,k>o,j+2k=D+i ^S- e j ak ~ *V 
Notice in particular the homogeneous case £(0; 0, 0) = eo — JC, C(0; y, 0) = — eo+i — JC. 



Proposition 3 For all y > 0, fixed q > and uniformly in s near s = 0, 

Cte v, q) = mjV~ lT ( s - \) < ( s ~ \ °> + ^ q y' 2s U^)+ 

+ — -y- s -* 22 a— r- K *-\ Vh + q) ■ 

r(a) tif^ \V^+^J 2 V y J 

§ Notice that the same consideration introduced previously about the variable q apply for the variable 
T, but recalling the discussion outlined in the introduction, we will assume positive definite temperature 
and will work out the T = limit once the right analytic continuation has been achieved. 
|| A further term ^— a appears at s = |. 
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Proof Let a k = X k + q. First, we isolate the (possible) vanishing terms: 

+00 

n=—oo k&K 

the second term needs no further comments, for what concerns the first, we proceed 
as follows. To start, assume y > to be fixed. In such a case, we provide a formula 
that can be analytically continued in s near s = 0, and we prove that such extension is 
uniform in y for all y > 0. We apply first the Mellin transform, 

+°o , ^00 +°° 

£ Y,l(yn) 2 + «J- = fh / E ^'E «-"<*. 

n=-oofceX ^ > J ° n=-oo keK 

and hence the multi dimensional Poisson summation formula, we get 



/•OO r> f— /-00 00 2 2 , 

■*o y \ ) Jo „_i i.cii' 



the integrals are known and give (where K u (z) is the Bessel function) 



' r f . - 1) C (. -i;o,0 + ^EE' (M ' 



«r(s) V 2/M 2 7 r(s) V\/afc 

The second term is a integral function of s, and it is also easy to see (using classical 
estimates for the Bessel functions) that convergence of the series is uniform in y for 
bounded y > 0. This means that the second term can be analytically extended for all s 
smoothly in y > 0. The first term can have poles, but dependence on y and s are clearly 
distinct. Using the results in Proposition and the known expansion for the Gamma 
function, we find out that there exists a closed neighborhood of s = where there are 
no poles independently from y, for all y > 0. This gives the thesis. □ 
In the following result is the main decomposition of the partition function of the physical 
model at any given temperature in term of the geometric zeta function. 

Corollary 1 

log Z(T,q) = #ker(-A M + q) logT-log Jj' (l - e "1? v/W ) + 

-^Res K(s; 0, q), 8 = -0 - ^(1 - log2)Re Sl (((s; 0, q), s = -0 + 

-ilogpC(0;2vrT,g), 
smoothly in T > and for all fixed q > 0. 

Proof From Proposition El we just have to work out the first term. Recall that 
_^ = 8 + 1S 2 + ( 8 3^ and r (s - |) = -2v^F - 2^(2(1 - log 2) - 7] a + 0(s 2 ), 
and use Proposition to write 

1 r. A C -l , , , ^/ 2\ ^1 , r. , , 2\ 



C [s - -; 0, q = — + c + cis + O(^) = — + i? + cis + 0(s 2 ), 
\ 2 / s s 

when s — > 0, where Ri = Res, (C(s; 0, q), s = — |). Then, f^yT (s — |) C (s — §; 0, q) 
-2y/nRi - 2^n[R + 2(1 - log2)i2i]a + 0(s 2 ), and the corollary follows. □ 



3. Low and high temperature expansions 

We state our results in the most general setting, using the zeta function introduced in the 
previous section. The expansions for the partition function at low and high temperature 
can be immediately obtained as corollaries (For an overview on the physical literature 
on the thermodynamic on the zero mode and the high and low temperature expansions 
the interested reader can have a look at ^1] and to the references there). 

Proposition 4 For y — > + , with fixed q > and uniformly in s near s = 0, 



C(s;y,q) 



-r (s - 1 V (s - i; 0, g) + 2IC q y~ 2 X R (2s) + O (e"*) . 



r( a ) 

Proof This follows immediately from Proposition El □ 
Corollary 2 For small T , and all fixed q > 0, 

log Z(T, q) = #ker(-A M + q) logT - ^Res (((s; -A M , BC), ~ ) + 
-(l-lo g 2)Re Sl (((s;-A M ,BC),~) ±+ 

V j,k>0,j+2k=D+l / 

Proposition 5 For y — > +oo, with fixed q > 0, and uniformly in s near s = 0, 
C(s; y, q) = ((s; 0, q) + T ( £±1 - s) ( R (D + 1 - 2s)y D ^ + o(y D ^ 



Proof Let = + q. We first isolate the (possible) vanishing terms as follows: 

oo 

C(s; M) = 2^E [M 2 + a k]- s + ((s- 0, g); 

n=l fceX 

let o- = ^, then: C (1) (s; 2/; <?) = y _2s 2 £n=i Efc 6 /J™ 2 + aa fc]~ s = S/~ 2 'Cr(s;?), and we 
study Co-(s;g). We start assuming Re(s) > si > ^y^, a G [<7i, 1]. Since convergence is 
absolute we can exchange summation indices as desired. Applying the Mellin transform 

1 W n=l fcei^ 

splitting the integral at t = -, in the first integral, at < 1, and hence we can use the 
heat kernel expansion for — Am, while the second integral gives a regular function of s 
for all s. For the first term, we get (recall a = y) 
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notice that we can not use the Poisson formula here, since we are interested in the small 
a expansion. Now, for each j, 

jxi .1 / '\ 00 roc 

/ CT t^i-^e-^dt = T ( s - a + 3 -\ Cr(2s - 2a + j) - I f+^^e^dt. 

n=l J° ^ ' n=l ^\ 

The first term can just have simple poles at a — | and a — In fact the Gamma 
function can have a pole only at non positive integer argument, but for even negative 
integers the zeta function vanishes. Thus, for each j and all s ^ a — |, a — 

lim 2 V (° t s+ ^ a - l e~ n2t dt = 2F (s - a + ^) ( R (2s - 2a + j), 

where the analytic extension of the above function of s is intended for Re(s) < a + |. 
For small ex, we can write 

r (*) Jo ^ u u n-s) 

and this is a regular function of s for s closed to 0, independently from a, and hence 
gives the desired result. The second term with x = at, and using the Poisson formula 
becomes 



OO 



a' 3 / x s - l 2j2e~^J2 e ~ Xak(X)dx = °~ S 0(e~°)i 

^ 1 n=l k=0 

for all s, since the second factor is a regular function of s for all s. Collecting and 
recasting the correct functions we get the thesis, where we put in evidence the leading 
term and use the Riemann's functional equation to rewrite its coefficient. □ 

Corollary 3 For large T , 

log Z(T, ,) = (£±i) Up + l)T D + o(T°). 

4. Finite temperature results for some particular geometries 

We turn now our attention to the case of a scalar field of rest mass q in two fixed 
geometries: a cubic box E D of edge I and a D-dimensional torus T D = S^ n x . . . x S^ n , 
at any fixed value of the temperature T. The operator A in the action is the negative 
of the .D-dimensional Laplacian A plus a constant potential in the Euclidean space 
and the difference in the two models is in the boundary conditions: Dirichlet boundary 
conditions on the boundary of the cube and periodic boundary condition for the closed 
domain. A complete system of eigenvalues is then: \ n ^(T, q, I) = (27rTn) 2 + ^-|/c| 2 + q, 
where n G Z is an integer and k = (ki, . . . , k D ) is a positive integral vector in (N ) D in 
the first case, but is any integral vector in Z D in the second one (21 is the length of the 
circles in the torus in the second case). The zeta functions are 



Ce D (s;27rT,g,/) = 



{27 rTnf + ^\k\ 2 + q 
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Cr D (s;2TrT,q,l) = ^ 



(n,fc)eZ D + 1 



71 



{2TiTnY + — \k\ 2 + q 



for Re(s) > -^p. Using the results of Section we just have to study the associated 
geometric zeta functions, namely 



Cb d {s; 0,q,l) 



fce(N ) D 



7T 



pl*l 2 + « 



Cr D (s;0,g,/) = ^ 



k& D 



71 



pW 2 + ? 



for Re(s) > y. This is the aim of this section, that is subdivided in three parts. In 
the first we introduce a multidimensional generalization of the Riemann zeta function 
useful to describe the geometric zeta functions, subsequently we calculate the main zeta 
invariants, and in the last part we give the partition function and the thermodynamic 
functions. From now on we will assume D > 1 if not otherwise stated. 



4-1. Multidimensional quadratic zeta functions 

Let q 2 be a complex constant that is not real and negative, n a D-dimensional vector 
with integer components in 7* D , and A a real symmetric matrix of rank D with positive 
definite associated quadratic form, then we can introduce the functions: 

ne(N ) D 

( D (s; Aq)=J2 ^ An + = <T 2s + n T An, q) 

nei D 

when Re(s) > y. Notice that the definition introduced for Cd(s; A, q) is ad hoc to avoid 
problems for the homogeneous case q = 0. In fact, while the definition given for £/j>(s; q) 
extends to q = 0, in the Eisenstein sum one must omit the null vector when q = 0. With 
the definition above, this simply means to omit the q~ 2s term in the homogeneous case. 
We will use the notation Cd( s ;<?) f° r CD(s;I,q) in the following. Notice in particular 
that £ (s; q) = q~ 2s and 2£i(s; 0) = Ci(s; 0) = 2( R (2s). Actually, the results we are going 
to give for the Epstein zeta functions hold true for a large class of zeta functions, that 
we introduce now. Let n G Z D and A be as above, b and x = (xj) be real D dimensional 
vectors, and assume < Xi < 1. Then, we define the functions: 

( D (s; A, b, x, q) = ^ K n + x f A ( n + x ) + hT n + <f]~ S , 
Cd(s; A, b, x, q) = ^ \( n + x) T A{n + x) + b T n + q 2 ]~ s , 

when Re(s) > y. Notice that q 2 must be non vanishing in the definition of 
Cd(s; A, b, x, q). A lot is known about these multidimensional zeta functions, in 
particular the homogeneous case, namely the Epstein zeta function, has been deeply 
investigated (see [HI] for a good overview, or locally cited references). Here, we collect 
a series of results that seem more interesting and useful for the present purposes. As we 
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will see, it is easier to get more general results for the Epstein series, and this is essentially 
due to the possibility of more effective use of the Poisson summation formula. All the 
proofs are based on classical techniques, namely the Mellin transform and the Poisson 
summation formula; since these tools were used for all the proofs in the previous sections 
we omit to give details here and refer the interested reader to the literature available 
for a deeper account on this subject. 

We begin by introducing some analytic representations. These will be useful to get 
all information about the analytic extensions of the zeta functions, as well as when 
calculations are involved to evaluate them at some particular value. 

Lemma 2 £ D (s; q) = -§&>-i(s; q) + } T (s - \) £d-i (s - ±; q) + 



V ' (n,/c)e(N ) D V / 



2^ 



For the Epstein type functions, we have the following lemma when q ^ 0: 
Lemma 3 



, / a „ n tt" r [ s 2) _o,i D 2n s n_ s ( Vn T A- 1 n\ - .— — ■ — — 



D 



- Y - f- - K D (ttVV - b T A-^n T A'^n), if 4q 2 - b T A^b > 0. 

VtetAr(s) n ^\ q 2 -\bTA-ib) s -^ v v 

When q = 0, we need some more notation. Let Ojj be the elements of A. Let the 
minor of a lt i in A, and Ai denotes the D — 1-column vector whose elements are the 
elements of the first line of Ai^. Let b\ and x\ be the D — 1 vectors whose elements are 
the last D — l elements of b and x, respectively. Let B be the D — l square matrix whose 
elements are bi-ij-i = dij — -^-^ , where the indices % and j run from 2 to D — 1. 



Lemma 4 < D ( S; a, o, *, o) = Ci(«; 01,1,0, *i,o) + ^t^Co-i {s - §; b, 0, x u 0) + j 



(«-iL 



X COS27T 



nixi H Ai(n + xi)ni 



(n T Bn)i( s -i) S 2 

C D (s; A, b, 0, 0) = Ci(s; ai,i, fti, 0, 0) + ^r(!) } to-i (« - ^ B > °> " 4^7 J + ri7) a M 
x E~=i E nez? -! cos (2Am + 6l ) ni] V )_ 



'™ TA l,l"-4S7T( 2A l™+ 6 l) 2 ) 5 " 



(^f^ n i+^- A i A i\ 6 ^-i h i A i,i 6 ^ 



where the last two representations hold if — ^-7 or — jbj A^bi are not negative integers. 
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More Chowla-Selberg type formulas [U] |3] can be found in [3T]- We also recall the 
important reflection formula [23 [HI] 

Lemma 5 7r- s T( S )( D ( S ; A, 0, 0, 0) = ^=Lv (§ - s) ( D (f - s; A~\ 0, 0, 0) . 

f -E^r (-f ) q D D odd 
Corollary 4 &(0; A, 0, 0, g) = { / f ' + 

I vfe^ Deven 
+ V (n T A~ 1 n)-Tir_i,(27rgv / n T A- 1 n). 

See also and [HI] for the Kronecker limit formula. 

Corollary 5 For n = 1, 2, 3, . . .; Co(-n; A 0, 0, 0) = 0; for n = 0, 1, 2, . . q ^ 

D odd 

Cij(-n; A, 0, 0, g) = q 2n + ( D {-n; A, 0, Q,q)={ (-i)-f! w # o n+D 



(n+f)l v / dHtA 



g,2n+ij ^ even; 



Eventually, using the representation introduced in the previous lemmas, or using 
classical methods, we get all information about poles, residua and particular values. 



Lemma 6 The function £r>(s;g) extends analytically to the whole complex plane up 
to simple poles at s = ~, . . . , ~ and s = — ~ — j, j = 0, 1, 2, . . wift residua: 

Re Sl = ^%E^E^ , 2W ^ ( ^ ) * = -A-s + 

(-1)" j_ V-f f D \ ■-! ■■'--/-• 



1, . . . , * £ 2N; £d(0; g) = i^- + i^- 



2j / r ^' +1 ) ' 



Notice that in the homogeneous case there are poles at s = ~, . . . , ^, with residuum 



Re Sl 0), s = f ) = f ^ j Tri , z = 1, 2, . . . , D, £ D (0; 0) 



2 D 



Lemma 7 Tne analytic continuation of Cd(s; A, 0,0, g) and Cd(s; A, 0, 0, g) are regular 
on the whole complex s-plane up to a simple poles at s = — — j, j ; = 0, 1, 2, . . . i/ £) 
odd and s = -j — j, j = 0, 1, 2, . . . , -j — 1 i/ D is even, respectively. The residua are, for 

both the functions, Resi (Cd(s; A, 0, 0, g), s = f - /) = ^^^^J— )' as &e / ore - 

f L> odd 

Moreover: Cd(0; A, 0, 0, g) = 1 + C(0; A, 0, 0, g) = \ { _^ qD 



— .D even 

T' 



Notice that in the homogeneous case the unique pole is at s = — with residuum 



jt, and C(0; A, 0, 0, 0) = -1 for all D. 



VdctAr. 
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4-2. Dirichlet boundary conditions 

In this case, the spectrum of —Am is positive definite, and we can write the geometric 
zeta function using the function q) just introduced: Cb d (s; 0, q, I) = |^£d (s; 
Notice that this extends continuously to the homogeneous case q = 0. Using 
Propositions El and El in Section and Lemma |H1 of the previous part, we get: 

Proposition 6 The function (B D {s',y,q,l) extends analytically to a regular function on 
the whole complex s-plane up to simple poles for all the s = P+ 2 1 ~ 3 , with j = 0, 1, 2, . . ., 
that are not non positive integers, with residua (k = —D — 1, —D, —D + 1, . . . , k (jt N): 

Re Sl [Cb d (s; y,q,l),s = --J = ^ - 




^ Ail 




i=0 j=0,2j-i-l=fc 

Notice that in the homogeneous case the poles are at s = ^j^, ^, ■ ■ ■ \ - 
Proposition 7 For all y > 0, fixed q > and I > 0, and uniformly in s near s = 0, 

Cb d (s; y, q, I) = ^- lW) v T ^--)t D {s--;—) + 

An 3 ,1 
v ' (n,fc)e(N ) I3 + 1 

As a corollary, we give the formula for the derivative at s — 0. Some care is necessary 
to deal with the first term when q ^ 0. In this case in fact, the function £d(s; q) has a 
pole at s — — |. Proceeding as in the proof of the corollary to Proposition El we get 

Corollary 6 For all fixed q > and I > 0, 

, , n n 27r2 n 1 V^V D \ (-1)* 

^(0;y,g,0 = ^(g,0--2to g J] (i-e"^), 

^ fce(N )o \ ' 

where A D (q, I) = -f [#0 + (2 - log4vr + Iog0i2i], Ri = Res, (s; ^) , * = -5) ■ 

Notice that in the homogeneous case the value at s = is for all D. Also, when 
q — 0, £d(s;0) is regular at s — — |; thus, an explicit formula for the constant when 
q = is Ad(0, Z) = — j^d (— 1;0), but it is more complicate otherwise. Using Lemma 
El if q = 0, when D = 1 we get twice the Riemann zeta function at s = —1, while 
for higher D only numerical evaluations are possible and give: £1 (— 1;0) = 1) = 
-i, £ 2 (-|;0) = 0.026127, £ 3 (-|;0) = -0.010015. Notice that these results can be 
easily obtained using the representations introduced in Lemma 121 since the series there 
converges very fast. The same computation can be made using the reflection formula 
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for the multidimensional zeta function given in Lemma but in that case the series 
converges very slowly and a much longer computation is necessary. 
We conclude this part with some remarks on the well know cases D = and D = 1. 
When D = 



-f Y^=2 [ 2 ) Cr(2j — In particular, when q = 0, the first term of the s-expansion 



z x {s; y, q) = ( Bo (s; y, q, vr) = ^(j/V + q)- s = q~ s + 2y~ 2s ^ (s; ^) ; 

for g ^ 0, zi(s;y,q) has a simple pole at s = \ — j, j — 0,1,..., with residuum 
r(ffy Zl(0;2/ ' g) = °' ^°'^) = "21og2sh^. For q = 0, ^(s^O)) = 
2?/~ 2s Cr(2s), has a single pole at s — ^ with residuum -, and zi(0;y,0) = — 1, 
z / (0;y,0) = 21ogy-21og27r. 

When D — 1, we have the zeta function associate to the Laplacian plus a constant 
potential on a cylinder [2~H] . 

oo 

Z2(s;y,q) = ( Bl (s;y,q,7r) = ^2^2(y 2 n 2 + k 2 + q)' 8 , 

neZ k=l 

and ^(Ojy.g) = -g, 4(0; 2/, g) = f J - (7- log2)f - 2 log ^=1 (l - e'^V 5 ^ - 
/ I 

near s = of such zeta function is well known 

4(0; y,0) = |i -21ogIJ (l -e->) = -2\og V (£j , 

where rj(z) is the Dedekind eta function. In such case, it is very easy to pass from the 
small y to the large y expansion using the well known modular transformation of the 
eta function, namely r\ (— ^) = \f^f]{ T )- It is also easy to see how the presence of a 
non homogeneous term breaks this symmetry. A deeper investigation about this point 
is performed in the next part. 

A final observation is about the large y expansion. As just noticed, this can be 
immediately worked out when D = 1. For higher values of D, it is harder. In fact, 
we can not use the analytic representation given in Proposition [7| for large y neither 
we know how the zeta function behaves under the modular transformation y — > -. 

a y 

Despite that, a direct approach is still possible for fixed D, consisting in expressing the 
zeta function in dimension D recursively in terms of the zeta function in dimension 
D — 1, and using the known behavior in dimension D = 1. Calculations are tedious but 
straightforward; the leading term in the expansion is consistent with the one obtained 
using Proposition but now we can get further terms. 

4-3. Periodic boundary conditions 

Recall that in the present case the spectrum of — A M is not always positive definite; 
more precisely, fC q — 1 if q = 0, but it vanishes if q > 0. We need the functions Cd and 
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(d of l4.1l in order to write the geometric zeta function: (rb( s > 0, q, Z) = \s\ 

q~ 2s + ^hCd (yS\ ^ j . Due to the presence of the first term, it is clear how this expression 
does not extend to the homogeneous case. In the present situation it is easier to deal 
with the two cases independently. This does not affect the poles, hence we can state the 
following unique result using Propositions |2 and El in Section EJ and Lemma l4l of 14.11 

Proposition 8 The functions (r D (s;y,q,l) and (r D (s;y,q,l) extend analytically to a 
regular function on the whole complex s-plane up to simple poles for all the s = —j, 
with j = 0, 1, 2, . . if D is even, and j = 0, 1, 2, . . . , : ^±1 — 1 ; if D is odd. The residua 

are for both the functions: Resi (Cr D 0; y,q,l),s= ^ - j) = .,Jp+i a S=tz- 

3 \ 2 3) tt 2 y 

Notice that in the homogeneous case, the unique pole is at s = with residuum 
D _ 1 ~ D+i — Ay' ^ n arLa lytic representation analogous to the one stated in Proposition 

0of the previous part is 



Proposition 9 For all y > 0, I > 0, and uniformly in s near s = 0, 

1 

2' vr 



Cr D {s;y,q,l) = " _, ^,'\ y L T [ s - ^ ) ( D [ s - ) + 



i 23 - 1 A - lr f iv ( 




47r s ^ > x * / n \ 2im /vr 2 IJIO 

K s-i —\-»\ k \+Q 



uniformly in q for q in any closed subset of the positive real axis, while 

i 23 - 1 / i\ * ( i 



Crz,(a;y,0,0 = ^-ir(> " lr V - 2 j Cd V ~ ^°)+ 2 y~ 2s Cn( 2s ) + 



T(s) V 2 ^ ^ \\k\J " s "t V ly 



™ =1 feezf 



Corollary 7 Uniformly in q > qo > 0: 

2tt 2 1 
(T D (0;y,q,l) = — 



£) even 

-D+l _ 

(-,)-a-i« i 



/ y ^ Vi 1 £odd, 



27T / 

<r D (0;j/>?>0 = ^(?>0— -2 log n (l 

Cr D (0; y, 0, Z) = ~R 1 - + 2( R (0) = -1, 
' 2/ 

(0; y, 0, Z) = £ D (0, Z) — + 2 log |- - 2 log J] 



e ! f 1 1 



w/iere Z) = -f [# + (2 - log4vr + logZ)iH R l = Res, (Cd (s; ^f) 

and similarly for the hatted ones. 



' S 2 1- 
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We get simple expressions for the constants when D is even or when q = 0. In particular, 
in the second case, B D (0, 1) = —j(d (— 0). Using LemmaEJand numerical evaluations: 
6 (-|;0) = 2Cb(-1) = -I 6 (-|;0) = -0.2286, ( 3 (-|;0) = -0.26493. 
Of particular interest is the case D — 1. As stated in the previous part, this is related 
with the Dedekind eta function, rj(z). Namely, assuming I = n for simplicity, 



lm \ [Cr 1 (°;2/,9, 7r ) + lo g<?l = 2 logy - 21og2vr - 41og?7 ( - 
This suggest to define the function 



V ( T , q) = -e mTB ^~> (l - e 2mrq ) JJ (l - e 2 ™ T V^+^j , 

n=l 

for real positive q and complex r with positive imaginary part. It is easy to check that 
lim g ^ + = ^ 2 ( r )- ^ * s a ^ so eas Y to realize that the presence of the non homogeneous 
term q, breaks modularity. On the other side, the modular transformation for the 
Dedekind eta function can be deduced using the symmetry in the definition of the zeta 
function (^{s', y, q 2 , 7r) under the exchange of the summation indices in the first term of 
the s-expansion near s = 0. Using the same symmetry for the function r)(iy, q), we get 
instead of the modular transformation the following relation 

(" \ 00 / 1 \ 2-7T °° / 1 \ 

-, 93/ J = log V (iy,q)-nq 2 y log y + 2-ny^ I | ) - 1),* + * £ ( f J Cfl(2j - l)(g?/) 2j . 

Since the behavior of r/(iy, g) for large ?/ is clear, the above expression can be used 
(exactly as it was for the Dedekind zeta function) to deduce the behavior for small y. 
We get, for y — > + : 



V(iy, l) = - — + -*q 2 y log y + 7r(g + 1) 
6y 



30 / I \ 



3=2 



2j i 2 
J -J- 7rq 



^Cfl(3) + 0(j/ 4 ). 



Thermodynamic functions 

We write now explicit formulas for the partition function of the models introduced in 
14.11 Such formulas can be used to get explicit expressions for all the thermodynamic 
functions. In particular, the behavior for low and high temperature are given. The 
partition function for a massive scalar thermal radiation at temperature T in a box of 
volume l D and on the torus To are, for any fixed positive I, 

logZ BD (T,q,l) = ±A D (q,l)±-log J] (l - e'Hf^A + 
i'^V D \ (-lY l 2i+1 q i+1 



logZ TD (T,q,l) = ±B D (q,l)±-log J] (l - e'^^A 



+ 



16 

these are smooth functions of the temperature T, for bounded T > 0, uniformly in the 
mass term q, for q in any closed interval of the positive real axis. When q = 0, 

iogz BB (T ) o,0 = — eD(~;o)i-iog n (i-e-^ |fc| )> 

V ' fce(N ) D 
log Z Td (T, 0, = {^y o) ± + logT - log JJ (l - e~^l) + I log p, 



where some values for the multidimensional Riemann zeta functions are given in 
14.21 and 14.31 respectively. Using these expressions, we get the behavior of the main 
thermodynamic functions: partition function, energy, entropy, pressure of the radiation, 
and specific heat. For low T, fixed q and I, we get on the box 

log Z Bo (T,q, I) = ^A D (q,l)-+ 

i-i) D+ \ i SV d \ (-i)M 2i +y +i L , 



D-ll 



S ED (r, 5 ,/)=0(r- 1 e-*), c BD (T, g ,() = 0(T- 2 e-*). 
while on the torus, we must distinguish the q = case: 

1 11 i f D even 

logZr D (T,q,l) = -B D (q,l)---logp-^ ^ ^ + 0(e"*), 

" '-1 u — 1 — 



, , I -D even 

= ~2fc(ff,/) + -logp<j Dodd +0(Te-*), 



Sr D (T,q,l) = O(T-V^), cr D (T,g,/) = 0(T"V 



T 



log z Td (t, o, = -^c o) i + io g r + i logp + o(e"*; 

F Td (T, 0, /) = -TlogT + C 0^ - ^Tlogp + 0(Te"^ 

S Td (T, 0, = logT + ~ logp + 1 + O^e-T), 
c TD (T,0,/) = l + O(T- 2 e-^). 
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Notice that the null mass case needs no independent treatment on the box. Recalling 
the remark at the end of 14.21 or the Proposition we get the behaviors for high T, fixed 
q and I (cfr 16 or 23J), 

log Z Bo (T,q,l) = -4^r f^±T) Cr(D+ l)T D + 0(T D ~^ 
log Z TD (T,q, I) = Ck(D + l)T D + 0{T D - 1 ). 

4-5. Critical volume 

As anticipate in the introduction, we will analyze in this section the dependence on the 
volume of the pressure of the radiation at finite temperature. The analysis is performed 
for the two models described in section 4, that only differ for the boundary conditions: 
periodic or of Dirichlet type. For simplicity, just consider the zero mass case. By 
definition, 

P(T,V) = -^T\ogZ(T,q,l) = lA( TC '(0;27rT,0,0) - ~(logp TC(0; 2ttT, 0, /)). 

Notice that applying corollary El for the Dirichlet boundary conditions, we always 
have the vanishing of the zero mass zeta function at s = 0, namely £(0; 27rT, 0, I) = 
for the box. For periodic boundary condition, we get a non trivial term involving the 
renormalization constant p. Thus the analysis in the following holds for the torus only 
if we assume the renormalization constant to be volume independent With this 
assumption, we get at each fixed temperature 

2vt t/ _d+i I \k\ 



^fcezf erv-b - 1 



\k\ 



D 



where the constants are 



H D = - l -C D [-l ), A- D = -i&(-i;0 

We show that for the periodic boundary condition there is, for all D, a value Vq 
of the volume where the pressure changes sign, being attractive when V < Vq. We 
also show that the same happens for the Dirichlet boundary condition at the physical 
dimension D = 3. For the torus, consider the function 



9d{x) = = + 9o(x), 

k& o e - - 1 e - - 1 

% As observed in the introduction, we are not going to analyze here the renormalization aspects of 
the model. Beside, notice that even with a renormalization parameter depending on the volume, the 
volume effect still exists, but critical volume depends also on the explicit form of p. This one, as pointed 
out in |14j . gives the connection between the model and the physical reality. 
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for x > and the following inequality 
Lemma 8 If x > ab and y > b (x, y, a, b > 0), then 
e ^y -A< (fy _ i) . 

Taking 6=1 and a = ^, we get the following bounds for gD^x), when x < xq, 



D i \ v D 1 ^ / \ 

< 9d{x) < - 



ex— l e * — 1 e 1+x o 



where Cd(xq) = JZha^D \k\e 1 + as o'*' is a positive constant. It is thus clear that the 

Sti, {0,l} 

pressure changes sign for some value of V (and fixed T) if the constant Hp is positive. 
To show that this is the case, just use the reflection formula (lemma 5) that for the 
function Cd(s; 0) takes the simpler form 

tt-T( s )Cd(s; 0) = tx s ~^Y - s] Co (^r-s;0 



.2 / V 2 

Using lemma |H] with 6 = 2 for the box, we get for the function 

Id{x) = = ^ + fo(x), 

fce (N )D e * - 1 e - - 1 

(x > 0) similar bounds when x < 

Vd Vd i _ . . 

< 9d(x) < + — — L D (x ), 

e x — 1 — 1 e 1+x o x — 1 

where Ld{xq) = J2ke(N ) D -{i} l^l e is a positive constant. Furthermore, using 

simple bounds for the norm, we get a bound for the constant Ljj(xq) 

iD / 2 \ D 

, e i+-o (2 

"2" . 



L D (xo) < 2 - / 2 n20 

1-e 



We can not prove that Kp is positive in general, but we can analyze explicitly the 
low dimensional cases. We find that Kp is negative for D = 2, but it is positive for D = 1 
and 3. This indicates that, up to renormalization, the presence/absence of a boundary 
does not affect the Casimir effect in the physical dimension for the model under study. 
Eventually, using the above bounds, we can provide bounds for the solution x = x* of 
/d(x) = Kb- For example, if D = 3, we get K 3 = 0.0050075, and with xq = 1.5923, 
C 3 (x ) < 0.23433, and 0.19684 < x* < 0.29613. 
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